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ABSTRACT 

We investigate how the conditions occurring in a protoplanetary disc may determine 
the final structure of a planetary system emerging from such a disc. We concentrate 
our attention on the dynamical interactions between disc and planets leading to or- 
bital migration, which in turn, in favourable circumstances, can drive planets into a 
mean-motion commensurability. We find that for a system containing a gas giant on 
the external orbit and a Super-Earth on the internal one, both embedded in a gaseous 
disc, the 2:1 resonance is a very likely configuration, so one can expect it as an outcome 
of the early phases of the planetary system formation. Our conclusion is based on an 
extensive computational survey in which we ask what are the disc properties (the sur- 
face density and the viscosity) for which the 2:1 commensurability may be attained. 
To answer this question we employ a full two-dimensional hydrodynamic treatment of 
the disc-planet interactions. In general terms, we can claim that the conditions which 
favour the 2:1 mean-motion resonance exist in the protoplanetary discs with mass ac- 
cretion rates in the range of 5-10^^ — 5 • 10^^ M0/year. For accretion rates higher than 
those needed for the 2:1 commensurability we observe a variety of behaviours, among 
them the passage to the 3:2 resonance, the scattering of the Super-Earth or the diver- 
gent migration caused by the outward migration of the gas giant. The results we have 
obtained from numerical simulations are compared with the predictions coining from 
the existing analytical expressions of the migration speed and the strength of the mean 
motion resonances. The conditions that we have found for the attainment of the 2:1 
commensurability are discussed in the framework of the properties of protoplanetary 
discs that are known from the observations. 
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1 INTRODUCTION 

Mean-motion resonances are ubiquitous phenomena in plan- 
etary systems. They carry important information about the 
formation processes and the further evolution of those sys- 
tems. Moreover, their occurrence may increase the chance 
of detecting planets which happen to be in a resonant con- 
figuration ( |Agol et al.(2005}[HoIman fc Murray|2005[ ). They 
might also be able to tell us more about the properties of the 
protoplanetary discs in which planets are born. The increas- 
ing number of low-mass planets being discovered in multi- 
pie planetary systems (|e.g. Mayor et al.|2009a (HP 40307) 



Mayor et al.|2009b (GJ 581)||Lovis et al.|2011 (HP 10180) 
is the first strong observational motivation for the studies 
of planetary formation and evolution. Another motivation 
is connected with the growing amount of data coming from 
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high angular resolution observations of protoplanetary discs 



(e.g. Andrews et al. 2009 Hughes et al. 2009 Isella et al 



2010 Olofsson et al. 2011 Huelamo et al. 20111. The ex 



citement in this area of research is stimulated also by the 
fast progress in theoretical investigations of the orbital mi- 



gration due to the action of tidal torques ( Lin & Papaloizou 



T993| |Ward|[T997l [Tanaka et al.|[2002l [Paardekooper fc Pa- 



paloizou 2008, 20091. The convergent migration of planets 



or protoplanetary cores is one of the most promising pro- 
cesses to explain the formation of resonant configurations. 
Extensive studies of possible resonant structures have been 



performed among others by|Nelson fc Papaloizou ( 2002 1 and 



Papaloizou fc Szuszkiewiczl ( 2005 20101. In Nelson fc Pa- 



paloizou ( 2002 1 the possible commensurabilities to be ex- 



pected in the system of planets in the Jovian mass range 



are investigated, while in Papaloizou fc Szuszkiewicz ( 2005 



20101 migrating planets in the terrestrial mass range are 



considered and the conditions for the occurrence of first or- 
der mean-motion resonances are provided. [Fogg fc Nelson| 
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Figure 1. The parameter space of the properties of the protoplanetary discs: surface density - viscosity plane. The central medium gray part 
denotes the region where planets attain the 2:1 conimensurability and the small hatched area favours the formation of the 3:2 resonance. 
Upper dark gray and lower light gray colours show the parameters for which the migration is divergent, and the white colour denotes the 
transition region where the variety of behaviours have been observed. 



(|2007L|Zhou et al.| ( |2005[ ), [Raymond et al.| ( |2006[ ), |Mandell 
|et al.| ( |2007[ ) have shown that terrestrial planets can grow 
and be retained in the hot-Jupiter systems and that their 
resonant captures by a gas giant are quite common. In the 
present paper we consider a system containing a Jupiter-like 
planet evolving together with a Super-Earth in a gaseous 
disc. We perform a survey of the difTerent disc properties in 
which this system is embedded using full hydrodynamical 
calculations to simulate the dynamical interactions between 
discs and planets. Our purpose here is to find out what are 
the conditions for the occurrence of the 2:1 conimensurabil- 
ity in such systems and at the same time to get an insight 
into the environment suitable for the formation of planets on 
the basis of the existence of certain mean-motion commen- 
surabilities in young planetary systems. It has been shown in 
Podlewska & Szuszkiewicz ( 2008 , 2009 ) that in the gas giant 



Super-Earth systems the planets embedded in a gaseous 
disc are easily locked in the commcnsurabilities if and only 
if the Super-Earth is located inside the orbit of the gas gi- 
ant. For this reason, we consider here exactly that kind of 
configurations in which the Super-Earth is the inner and the 
gas giant is the outer planet of the system. In the previous 
paper (Podlewska & Szuszkiewicz 20081, where the same 



system has been studied, the 2:1 commensurability was not 
included because of computational limitations. 

This paper is organized as follows: In Section [2] we de- 
scribe how our computational survey has been performed. 
Next, we illustrate the relevant features of the evolution of 
Jupiter-like planets in discs with different viscosities. In Sec- 
tion H] we give the results of our survey presenting differ- 
ent configurations achieved during the tidal evolution of the 
planets in the disc. These results are compared in Section 
[5] with the theoretical expectation based on the analytical 
and semi-analytical expressions for the migration speed of 



the planets and the strength of the mean motion resonances. 
In Section [6] we point those observed protoplanetary discs 
in which the physical conditions that favour the occurrence 
of the 2:1 resonance are present. Finally, our findings are 
summarized and discussed in Section [7] 



2 SETTING THE SCENE FOR THE SURVEY 

In order to determine possible resonances induced by planet 
migration in the presence of a gaseous disc in a system con- 
taining a Super-Earth and a gas giant, we have performed 
a computational survey by means of full two-dimensional 
hydrodynamic simulations. The mass of the Super-Earth is 
taken to be 5.5 M^ and that of the gas giant 1 Mj (here 
M0 and Mj denote the masses of the Earth and the Jupiter 
respectively) . This pair of planets is used in our simulations 
as a probe of the disc properties which favour the resonant 
capture of the Super-Earth by the gas giant. We have fixed 
the disc aspect ratio to h — 0.05, which is a typical value 
for protoplanetary discs. Herewith, we have assumed for the 
initial surface density the classical power law parameteriza- 
tion 



E(r) 



^" {5:^) 



(1) 



with p — 0. We will comment on these assumptions in Sec- 
tions [5] and [6] In our survey different discs will be analyzed, 
whose properties will be distinguished by using two param- 
eters, namely the surface density value E and the coefficient 
of kinematic viscosity v, which is constant in space and time. 
In the following, these two parameters will be referred to as 
"disc properties" . The parameter space for our studies is il- 
lustrated in Fig. [1] The value of the surface density E varies 
from 10"* to 10~^ in the dimensionless units used in the 
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code. The unit of mass is the mass of the central star M 
(here we take M = Mq) and the unit of length Vp is the 
initial position of the inner planet (the Super-Earth in our 
case). The unit of time in the code is given by {GM/rpY'^ 
(G is the gravitational constant). This quantity amounts to 
(l/27r) times the orbital period of the initial orbit of the 
Super-Earth. The value of 6 ■ 10"'* of the surface density 
in these units is consistent with the Minimum Mass Solar 
Nebula (MMSN) ( ,Hayashi„1981J at the present position of 
Jupiter in our planetary system, i.e. 5.2 AU. The range of the 
surface densities considered here expressed in physical units 
is the following: 30 - 300 g/cm^. In order to constrain the 
values of the kinematic viscosity v for our studies, we have 
used the relation between the kinematic viscosity and the 
accretion rate M observed in T Tauri stars, which are the 
best laboratories for studying the early stages of planetary 
evolution. For a stationary thin accretion disc this relation 
reads 



T.v = 



M 
3^' 



(2) 



The upper and lower curves in Fig. [T] are the constant ac- 
cretion rate levels. For the upper curve we have adopted 
3 • 10"'' Mo/year and for the lower one 3 • 10~^ M0/year, 
which are typical accretion rates observed in T Tauri stars. 
The viscosities obtained in this way vary from 6 ■ 10^'' to 
6 ■ lO"'^ in terms of the inverse of the Reynolds number at 
the initial location of the Super-Earth. Those two curves, 
together with the minimal and maximal values of the sur- 
face density, form the region of parameter space which is of 
interest for our studies. In the region defined in this way we 
have chosen 47 initial disc configurations for which we have 
performed numerical simulations. Each black dot in Fig. [l] 
represents a protoplanetary disc of given surface density and 
viscosity. The circles around some of the black dots indicate 
those discs which are of particular interest and are described 
in detail in the next Sections. The sampling of the parame- 
ter space is not uniform. We have calculated the evolution 
of planets only for those discs which were relevant for our 
purpose, namely for determining the variety of the final ar- 
chitectures of the planets. 

We have performed full two-dimensional hydrodynamic 
calculations using the grid based code NIRVANA. For details 
on the numerical scheme and the adopted code see [Nelson] 
et al. (20001. The planets interact with each other and with 



the protoplanetary disc in which they are embedded. These 
dynamical interactions take place in the central potential 
of a solar mass host star. The disc is locally isothermal. 
It extends from 0.33 to 4 in dimensionless units. We use 
the resolution of 384 x 512 grid cells in the radial and the 
azimuthal directions respectively. The grid spacing in both 
coordinate directions is uniform. The computational domain 
is a ring, in which the azimuthal coordinate 4> takes its values 
in the interval [0, 2ti\. We choose open boundary conditions 
in the radial direction and periodic in the azimuthal one. 
We do not exclude any material from the Hill sphere and we 
discuss the consequences of this in Section |5] The potential 
is softened with the parameter 6 = 0.8H, where H is the 
semi-thickness of the disc at the planet position. The initial 
eccentricity of both planets is set to zero. 
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Figure 2. The semi-major axis evolution of the migrating 
Jupiter-like planet in the disc with different values of the vis- 
cosity and with the surface density corresponding to the MMSN 
at 5.2 AU. The curves are labeled by the value of the viscosity in 
the disc given in terms of the inverse of the Reynolds number. The 
'viscous time' means the time expressed in viscous time units. 
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Figure 3. The surface density profiles at 0.05 t^ for three dif- 
ferent viscosities: 5- 10~® (grey line), 10~^ (black line) and 10"'* 
(dots). 



3 GAS GIANT MIGRATING IN DISCS WITH 
DIFFERENT VISCOSITIES 

It has been established that the migration of the gas giant in 
protoplanetary discs depends on the value of the viscosity 



( Lin & Papaloizou J1986, ,1993| ) . The viscosity controls the 
processes responsible for the gap opening in the disc which 
is crucial for further evolution of the orbit of the planet. 
We have checked that for viscosities lower than 10~^ and 
typical surface densities, the migration of the Jupiter-like 
planet in a reasonably good approximation, proceeds on a 
viscous timescale as in the classical Type II regime. 

Relatively recently, |Crida fc Morbidelli| ( |2007| ) found 
that for sufficiently high values of the viscosity, correspond- 
ing to Reynolds numbers of around 10'', the migration of 
the gas giant tends to change its direction and the planet 
moves outwards. The values of viscosity for which they have 
obtained the outward migration of the Jupiter-like planet 
are present also in our parameter space so, we have been 
able to verify if the gas giant behaves similarly in the discs 
analyzed in our survey. The results, illustrated in Fig. |2J 
show the evolution of the semi-major axis of the Jupiter 
mass planet placed in the protoplanetary disc with the sur- 
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face density equal to the MMSN at 5.2 AU and different 



kinematic viscosities. Following Crida & Morbidelli (20071, 



in Fig. [2] and in all other figures where time is present, the 
time is measured in units of the viscous time t^, which is 
defined in the following way 



t. = 



' p 

V 



(3) 



If the viscosity becomes higher than 10~^, the Jupiter mi- 
gration slows down, and if the viscosity reaches the value of 
5- 10~^ or higher, the migration of the Jupiter is reversed. It 
means that the Jupiter-like planet will not be subjected to 
proper type II migration for sufficiently high disc viscosities. 
Our results are thus consistent with those obtained by | Crida] 
& Morbidelli ( 2007 1 despite the differences in the initial sur- 



face density profile and the numerical treatment adopted in 
this paper. To complete the picture of a gas giant migrat- 
ing in the discs with different viscosities, we present also the 
surface density profiles (Fig. |3| for three of our simulations 
characterized by viscosity values of 5 • 10~^, 10~^ and 10^* 
respectively. All profiles are shown at the same time of the 
simulations, namely at 0.05 t^. This figure clearly illustrates 
how the gap profile depends on the value of the viscosity. It 
is also worth mentioning that the mass of the inner disc crit- 
ically depends on the position of the gas giant with respect 
to the inner edge of the grid as has been shown by |Crida| 
& Morbidelli ( 2007[ ) . The outward migration of the Jupiter 
has an important consequence for the resonant captures de- 
scribed in Section [4] 



4 EXPECTED COMMENSURABILITIES FOR 
A VARIETY OF PROTOPLANETARY DISC 
PROPERTIES 

The planets have been placed in various discs with surface 
densities and viscosities taken from the range shown in Fig. 
fTland described in Section [2] The initial separation between 
the planets has been chosen to be slightly larger than that 
required for the 2:1 resonance. Next, we have followed the 
evolution of this system for roughly 2 • 10^ orbits of the inner 
planet i.e. the Super-Earth. Depending on the final configu- 
ration of the planets, we have divided the whole parameter 
space into five regions. The medium grey part of the dia- 
gram in Fig. [1] is characterized by the fact that the most 
likely outcome of the evolution of the migrating planets is 
the resonant capture in the 2:1 commensurability. The disc 
with properties lying in the small central hatched region 
favours the formation of the 3:2 resonance. For discs with 
viscosity smaller than 4 • 10~® the differential migration is 
divergent and we do not expect to find any commensura- 
bility in this case. The same is true for discs with surface 
densities and viscosities taken from the dark grey part of 
the diagram in Fig. [T] In this region the migration of the 
Jupiter is directed outward, while the Super-Earth migrates 
inward, so that the resonant capture is not possible. A tran- 
sition region in which a variety of different behaviours has 
been observed is denoted in white. 

At this point, we would like to describe briefiy the typ- 
ical behaviour of the system in each region of the parameter 
space shown in Fig. IT] 



4.1 The 2:1 commensurability 

The capture in the 2:1 resonance is a very likely outcome of 
the planet evolution in the disc. Indeed, the medium dark 
grey region of the diagram in Fig. IT] in which the physical 
conditions in the disc favour the achievement of just that 
commensurability, is the most pronounced one. As it turns 
out from Fig. Jl] the best conditions for the occurrence of 
the 2:1 resonance exist in the region in which the surface 
densities are relatively low (from 10~* till 3.5-10"'') and the 
viscosity ranges from 4-10~^ till almost 10~^. For higher 
surface densities, in the range from 4.5-10"'' till 10""^, the 
spread of viscosities useful for attaining the 2:1 resonance 
shrinks to the interval from 4-10"'' till 10"^. 

First, we investigate how the planets locked in the 2:1 
resonance evolve in time in discs which differ among them- 
selves only by the value of the viscosity, while the surface 
density is kept fixed. For this purpose we choose E = 3 - 10"'* 
because, as it is shown in Fig. [T] this value of the surface 
density gives us the biggest range of viscosities for which the 
resonance under consideration is possible. From our analysis 
it turns out that the semi-major axis ratios oscillate always 
around the position of the resonance and we can identify two 
characteristic behaviours of this oscillation. The first kind of 
behaviour appears if the viscosity is 3- 10"'' or lower. In that 
case the semi-major axis ratios oscillate with a large ampli- 
tude, the same for all discs. The second kind of behaviour 
occurs if the viscosity is higher than 3 - 10"^ then the am- 
plitude of the oscillations decreases together with increasing 
viscosity. Examples of these two behaviours are displayed in 
Fig. p] (left panel) for two values of the kinematic viscos- 
ity, namely 5 - 10"^ (black) and 3 - 10"® (grey). The am- 
plitude of the oscillations is a straightforward consequence 
of the excitation of the Super-Earth eccentricity, shown in 
Fig. J4] (right panel). Let us also notice that after the capture 
in the mean motion commensurability the resonant interac- 
tions drive the eccentricities of the Super-Earth until their 
action is balanced by the circularization due to disc tides. 
For all discs with viscosities 3 - 10"® or lower we have ob- 
served that the eccentricities of the Super-Earth increase 
until they attain an equilibrium value of about 0.62. An ex- 
ample is given in Fig. J4] (right panel, grey line). For discs 
with very high viscosities there is still a lot of material in 
the vicinity of the Super-Earth, so the circularization forces 
act more efficiently and the final eccentricity (black curve in 
the right panel of Fig. B| grows to a value lower than 0.62. 
Thus, there is a correlation between the final eccentricity 
and the surface density of the disc around the position of 
the Super-Earth. In the case of higher surface density the 
eccentricity is lower. 

Next, we perform a similar investigation but this time 
we fix the value of the viscosity in the disc and check how the 
evolution of planets locked in the 2:1 resonance depends on 
the initial surface density of the disc (or equivalently on disc 
mass). The value of viscosity is taken to be i/ = 10"®, which 
is a typical value for protoplanetary discs. In Fig.j5]we show 
the evolution of the semi-major axis ratios and the eccentric- 
ities of the Super-Earth for discs with this typical viscosity 
and two values of the surface density, namely E = 10"'' 
(grey colour) and E = 6 - 10"'' (black colour). The value of 
the oscillation amplitude of the semi-major axis is again the 
consequence of the eccentricity evolution. The eccentricity 
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Figure 4. The evolution of the semi-major axis ratios (left panel), the Super-Earth eccentricities (right panel) due to the disc-planet 
interaction for two discs with the same surface density equal to 3 • 10"'* and two values of the kinematic viscosity, namely 5 ■ 10~^ (black) 
and 3 • 10~^ (grey). 
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Figure 5. The evolution of the semi- major axis ratios (left panel), the Super-Earth eccentricities (right panel) due to the disc-planet 
interaction for two discs with the same viscosity equal to 10~^ and two values of the surface density, namely 6 • 10~* (black) and 10"* 

(grey). 



of the Super-Earth in the 2:1 resonance is bigger if the ini- 
tial differential migration of the planets is faster which is in 
agreement with the Eq. (16) of Crida et al. (20081. It can 



be seen from Fig. [5] that both presented there eccentricities 
do not reach an equilibrium value. For the disc with higher 
surface density (E = 6 • 10~*) and viscosity v — 10~^ it 
occurs, because the Super-Earth arrives at the inner edge of 
the disc before the eccentricity attains an equilibrium value. 
Instead, for the disc with lower surface density (E = 10"*) 
and the same viscosity there is not enough material around 
the position of the Super-Earth in order to balance the res- 
onant interactions. This is why the eccentricity continues to 
grow till the end of the run. 

In all simulations in which the 2:1 resonance has been 
observed the angles between the apsidal lines and the reso- 
nant angles librate around zero. In Fig. [6] we plot the evolu- 
tion of the angle between the apsidal lines (black colour) and 
the resonant angle, defined as 2Aj — As —uJs (grey colour) for 
E = 10"* and v = 10~^. Here, \j and As are the mean lon- 
gitudes of the Jupiter and the Super-Earth respectively. The 
longitude of the pericentre of the Super-Earth is denoted by 
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Figure 6. The evolution of the angles between the apsidal lines 
(black colour) and the resonant angle (grey colour) for the disc 
with the kinematic viscosity 10^^ and the surface density equals 
to 10"^. 



there is no possibility of attaining the commensurability be- 
cause the differential migration becomes divergent. 



Our calculations have been able to determine also the 
lower limit of the viscosity necessary for the occurrence of 
the 2:1 resonance. We have noticed that below v = 4, ■ 1Q~^ 



4.2 Other configurations 

In the previous Subsection we have singled out within the 
parameter space of typical physical properties of protoplan- 
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Figure 7. The evolution of the semi-major axis ratios (upper left panel), the Super-Earth eccentricities (upper right panel), the angle 
between the apsidal lines (bottom left panel) and the resonant angle (bottom right panel) due to the disc-planet interaction for a discs 
with the surface density: 6 ■ 10"'* and the kinematic viscosity equal to 2 ■ 10~^. 



etary discs the region which is favourable for the achieve- 
ment of the 2:1 resonance. We have also characterized the 
final configurations of the planets in this commensurability. 
However, the planet evolution in gaseous discs does not al- 
ways end up in the 2:1 commensurability. In the following 
we would like to describe which are the alternatives to this 
behaviour and to determine the disc properties that lead to 
these different possibilities. 

Let us start from the region in which the planets pass 
through the 2:1 resonance and approach the 3:2 commensu- 
rability. This is the small hatched region of Fig. [l] in which 
the surface densities vary between 4.5 • 10"'' and 6 • 10~* in 
dimensionless units and the viscosities take their values in 
the narrow range 2-3- 10"^ For E = 6-10"* and z/ = 2-10"^ 
the planets enter the 3:2 commensurability. The semi-major 
axis ratio oscillates around the value 1.32 (Fig. [7|. The ec- 
centricity of the Super-Earth reaches a value of 0.45 which is 
lower than what we get for the 2:1 commensurabilities. This 
result is in good agreement with that previously presented in 



Podlewska & Szuszkiewicz (20081 for the 3:2 resonance. The 



angle Ao; between the apsidal lines and the resonant angle 
librate around the value of 230 degrees (4 radians). Usually 
one observes that the equilibrium value of Auu is either 
or 180 degrees. However, it has been shown in |Beauge "et| 
|al.| ( |2003[ ) that Aw is described by a function of the masses 
and the eccentricities of the two planets and that, for high 
eccentricities, it can librate around the values that are dif- 
ferent from or 180 degrees (e.g. |Kley et al.] (j2004)). In a 
disc with the same surface density (E = 6 • 10~ ) but with 
higher viscosity {u = 3 ■ 10~^), the 3:2 resonant trapping is 
also present, but it lasts only for roughly 470 orbits and after 
that the planets scatter. At the moment in which the planet 



scattering takes place, the eccentricity of the Super-Earth 
equals 0.38. The 3:2 commensurability has been achieved 
also in a disc with the surface density E = 4.5 ■ 10~* and the 
viscosity u = 3 ■ 10~'. In this case, the eccentricity of the 
Super-Earth is equal to 0.3. The system does not abandon 
the 3:2 commensurability till the end of our simulations. For 
the same surface density E = 4.5 • 10~* and lower viscosity 
{u = 2 ■ 10^'') we have again scattering in the system, but 
this time it occurs after staying for 400 orbits in the 3:2 
resonance. 

Below the threshold of the viscosity necessary for the 2:1 
resonant trapping, which occurs at 4 ■ lO"'^, the migration 
is divergent, see the light grey part in Fig. [l] This region 
of parameters does not contain resonances, so we will not 
discuss it any further. 

More interesting is the white region of high viscosities 
of Fig. [T] In this part of the diagram we observe unstable 
behaviours leading to major changes in the planet configu- 
rations. One of the examples of those instabilities is shown 
in Fig. [8] It is the case of a disc with surface density equal 
to 4.5 • 10"'' and viscosity 5 • 10~^. The gas giant evolving in 
the disc captures the Super-Earth in the 5:3 mean-motion 
resonance, as it is indicated by the fact that the relevant res- 
onant angle is librating around a constant value (see lower 
right panel of Fig. jsl. At a certain time after the instant 
0.1 in viscous time units the low-mass planet starts to go 
around its host star on a clearly chaotic orbit. Further evo- 
lution leads to the orbit crossing and the gas giant becomes 
the internal planet in the system. Moreover, the final config- 
uration after the orbit crossing is close to the 1:2 resonance. 
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Figure 8. An example of unstable second order comniensurability 5:3, which ends up with orbit crossing and evolves into a configuration 
close to the 1:2 resonance. The surface density is equal to 4.5 ■ 10"'* and the viscosity has the value of 5 ■ 10~^. The upper left panel 
shows the semi-major axis beheviour for both planets. The upper right panel illustrates the changes of the semi-major axis ratio in time. 
The lower left panel displays the evolution of the eccentricity of the Super-Earth orbit. Finally, the lower right panel presents the angle 
between the apsidal lines. 



5 THEORETICAL ASPECTS OF THE 
RESONANCE SURVEY 

In Fig. [l] we have shown the behaviour of the planets in the 
surface density-viscosity plane obtained from hydrodynami- 
cal simulations. It might be useful to compare our numerical 
results with those obtained from the analytical expressions 
for the migration speed provided in the literature. First, 
we assume that the Jupiter migrates according to the type 



II migration (Lin & Papaloizou (19931; Crida & Morbidelli 



(2007 1 ) and its migration time can be estimated by 



TII = 



2Tn 



(4) 



where rj denotes the radial position of the gas giant, ra,j is 
the mass of the gas giant, and Vu is the torque acting on 
the Jupiter. The torque is given by the following expression 
( Crida fc Morbidelh|2007| ) 

Vn = 27rr+Ez,p(r+)«4r+)\/GAfr+ (5) 



where 



T,Lp{r,t) = EoT 



-5/4 



■v^-V^' 



inf 



^T 



e^p{-a{GMrf /T) 



is a surface density profile taken from |Lynden-Bell fc Pringle| 
( 19741. Here, T = 12avt + l is a function of kinematic viscos- 
ity V and time i, Rin} is the inner edge of the disc, Eo and a 
are constants defining the initial shape of the surface density 
distribution. The radial velocity Vr as well as all other quan- 



tities in the expression for the torque acting on the planet 
is evaluated at the external edge of the gap opened by the 
gas giant, namely at r'^ — rj + 2Rh where Rh is the radius 
of the Hill sphere. Having specified the migration time for 
the gas giant we turn our attention to the Super-Earth. The 
evolution of the Super-Earth proceeds according to the type 



I migration as described in Tanaka et al. ( 2002 1 



Ti = (2.7+ l.lp)-^ — -^/I'n. 

TTls 2-isTs 



(6) 



where p denotes the surface density slope of the disc (see 
Eq. ([l|), JTis is the mass of the Super-Earth and Q,a is the 
angular velocity at the Super-Earth location r^. Thus, we 
can easily predict when the differential migration will be 
convergent, which is essential for the occurrence of the mean- 
motion resonance. From the Eqs. ([5| and ([6| the criterion 
for the convergent migration reads 

(1.364 + 0.541p) 
' > ^ 3^ ^" 



K^[l- ^{GMr+Y 1- 



1/2 



(7) 



Moreover, we can combine the expected differential migra- 
tion speed with the strength of the 2:1 and 3:2 commensu- 
rabilities given by |Quillen| ([2006|, [Mustill fc Wyatt| ( [20TT| 
in order to guess if the capture will take place or not. 
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The resonant capture for the first order resonances in 
the restricted three body problem occurs when 

' ^ '" (8) 






T]crit 



n 



where rjcru is the critical mean motion drift rate and fi is 
the angular velocity of the Jupiter. In the case of internal 2:1 
resonance ricrit ~ 22.7(mj/M)*'^ ar id for the 3:2 commen- 
surability ijc^t = 126A{mj/Mf^^ ( |Quillen||2006| ) . Finally, 
we make an attempt to draw the theoretical analog of our 
Fig. [I] The results are shown in Fig. [9] (left panel). In this 
figure we have used the same grey scale and patterns as in 
Fig. IT] in order to facilitate the comparison. 

The overall theoretical picture is similar to the one ob- 
tained in our numerical simulations. However, the details 
do not match exactly. The convergent migration takes place 
for viscosities higher than 1.2 • 10^^. This gives the region of 
the divergent migration smaller than in Fig. IT] The region of 
the 2:1 resonance occurrence in both numerical and theoret- 
ical approaches occupies almost the same place in the space 
of the disc parameters. The only significant differences are 
present for relatively low surface densities and high viscosi- 
ties. The hatched region with high viscosities (see Fig.lo] left 
panel) for which there are the conditions for the 3:2 com- 
mensurability, is much more extended than the analogous 
region that results from our numerical simulations. In the 
white region we expect the 4:3 commensurability. However, 
we would like to stress here that the picture discussed above 
has been obtained under the assumptions that the planets 
evolve according to classical type I and type II migrations. 

In a more realistic situation, when the gap opened 
by the gas giant is not deep and clean the torque acting 



on the planet needs to be modified ( Crida & Morbidelli 



(20071). Therefore, in a more refined approximation, we have 
changed the assumption about gas giant migration. Namely, 
instead of applying Eq.(lO) of [Crida fc Morbidelhl ( |2007| l as 
we did in order to obtain the diagram displayed in Fig. [9] 
(left panel) we use Eq.(15) of that work, so now the torque 
reads 



F = 27rr+ELp(r+)?;r(r+)\/GMr+ x 
WRh ^jr/ vr{rj) 



i-fiP) 



■^r+ VGMr+«r(r+) 



f{P) 



dlogS/5 
dloer 



(9) 
where E is the density inside the gap, B is the second Oort 
constant. The function 



f{P) 



(P-0.541)/4 ffP< 2.4646 

1 - exp(-P"-^V3) if P ^ 2.4646 



describes the gap depth expressed as a ratio of the gap sur- 
face density to the unperturbed density at r"*" and P is de- 
fined by 

iH 50 



P = 



ARt 



+ 



{mj/M)R 



<1 



where R is the Reynolds number. In this way we are able 
to take into account the torque exerted on the outer disc by 
the gas in the gap and the corotation torque. The result- 
ing picture can be viewed in Fig. p\ (right panel). As it is 
possible to see in this figure the hatched area has shrunk 
towards higher values of the surface density, which is in bet- 



ter agreement with the numerical findings. The region of the 
divergent migration for low viscosities remains unchanged. 
This is because, for low viscosities, the torque expressed by 
Eq. ([5| is very similar to that of Eq. ([9|, so the line divid- 
ing the regions of the divergent and convergent migrations 
can be obtained from Eq. ([7| also in this case. Moreover, 
the region of the divergent migration has appeared for high 
viscosities as expected. Whenever the migration is conver- 
gent, the planets will be captured either in the 2:1 or 3:2 
resonances. The line between the regions of the occurrence 
of these commensurabilities can be found using the criterion 
Eq. ([8| in which tu is expressed by Eq. Q, but Fu is re- 
placed by the torque given by Eq. H. The stability of the 
resonances requires further studies. Other methods must be 
used in order to determine the final outcome of the migra- 
tion of the resonant pair of planets. The numerical approach 
demonstrates to be very useful to tackle such problems and 
enables us to identify those regions in the parameter space 
in which the commensurabilities are not only form, but also 
are maintained at least for the time accessible to full hy- 
drodynamic calculations. In Fig. [l] those regions has been 
indicated by the medium grey colour and the hatched pat- 
tern. The rest of the parameter space, where the migration of 
planet is convergent has been left in white. The whole variety 
of behaviours has been observed in this white region includ- 
ing transient second order resonances and the orbit crossing, 
which indicates that not all resonant captures shown in Fig. 
[9] are stable. 

As already explained, in order to draw Fig. |9] we have 
used the assumption that the Super-Earth migrates ac- 
cording to [Tanaka et aL\ (|2002[), but recently it has been 
shown by Paardekooper & Mellema (20061, Baruteau & 



Masset (20081, Kley fc Crida (20081, Paardekooper et al 



( 2010 2011 1 that in non-locally isothermal discs the rate and 



even the direction of the type I migration can be strongly 
changed. Therefore, we discuss briefly how our results would 
be affected by considering non-locally isothermal discs. If the 
Super-Earth migrates inward slower than in our simulations 
we would observe a faster differential migration and in con- 
sequence the regions of occurrence of the 2:1 resonance in 
our surface density-viscosity plane would be shifted down- 
ward. In the situation, when the Super-Earth migrates in- 
ward faster than in our calculations, it would be more dif- 
ficult to obtain the convergent migration. Thus, the regions 
of occurrence of the 2:1 resonance would be shifted upward. 
However, we expect that the global picture should be simi- 
lar. In the case of the Super-Earth migrating outwards, the 
most likely effect would be the scattering of the low-mass 
planet. 

We have also checked the influence of the surface density 
proflle on our final results. In our calculations we have used 
the initial density profile with power-law index p = 0, while 
the observations suggest 0.5 < p < 1.5 ( Pietu et al.||2007 |. 
In doing that, we have run simulations with p = 1 taking 
E = 6 ■ lO""* and viscosities u = 10"^ 10"^ 2 ■ 10"^ 3 • 
10~^ and 8 • 10~^. It has been found that the differential 
migration is slower than for p = 0, so that the global picture 
is qualitatively the same but the regions of the characteristic 
behaviours in Fig.fTlare moved upward. For example the case 
with E = 6- 10~* and v — 2- 10~® gives the 3:2 resonance for 
the disc with p = 0. In the same disc with power-law index 
p = 1 we have obtained the 2:1 commensurability. 
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Figure 9. The parameter space of the properties of the protoplanetary discs: surface density - viscosity plane obtained under the assump- 
tions that the planets evolve according to the classical type I and type II migration respectively (left panel) and that the evolution of the 
Super-Earth is the type I migration and that of the Jupiter is modified by the gas possibly present in the gap | |Crida &c Morbidelli|2007l l 
(right panel). 



Moreover, as we have mentioned in Section [2] we do 
not exclude any material from the Hill sphere. To justify 
our approach, we have performed numerical simulations in 
order to investigate the influence of the material close to the 
planet on the gas giant migration rate. We have observed 
that excluding the material from the Hill sphere results in 
the slowing down of the migration of the Jupiter, which has 



been already shown by Crida et al. (20091. Thus, the differ- 
ential migration is also slowed down and the global picture 
would be the same as taking higher p. That means that the 
region of the resonant behavior will be shifted toward higher 
viscosities in our parameter space. 



6 OBSERVATIONAL ASPECTS OF THE 
RESONANCE SURVEY 

In our calculations we have started the investigations by 
choosing the masses of the planets and setting the main 
properties of the protoplanetary disc in which they are em- 
bedded, namely its surface density distribution, aspect ratio 
and viscosity. Ideally, one would like to use observed quan- 
tities in order to fix the initial conditions for such calcula- 
tions. However, regardless of the fact that circumstellar discs 
are the natural outcome of the star-formation process, their 
structure is still not known in sufficient details in order to 
set such conditions using direct measurements. 

The most powerful technique which could be able to 
change this situation and provide the most common proper- 
ties of protoplanetary discs, is that of high resolution imag- 
ing. Recently, [Isella et al.1 ( |2010[ ) have presented the CARMA 
(Combined Array for Research in Millimeter-wave Astron- 
omy) observations of the thermal dust emission from the 
discs around RY Tau and DG Tau with the best resolution 
achieved to date at wavelengths of 1.3 mm and 2.8 mm, 
namely 0.15" (which gives 20 AU at the distance of the 
Taurus cloud) . On the basis of their high angular resolution 
observations, Isella et al. (20101 have determined the disc 



vary by almost an order of magnitude. In our work, as men- 
tioned in Section|2] we simply assume the classical power law 
parameterization of Eq. [l] with p = 0. A better assumption 
for the surface density distribution will be possible only with 
the future advances in observational techniques. Taking the 
values of the surface density derived by Isella et al. (2010) 
we have placed DG Tau and RY Tau in the diagram of Fig. 
[l] The result is displayed in Fig. [lO] which shows the same 
regions as Fig. IT] DG Tau is located in that region of the 
parameter space in which a 3:2 resonance might be achieved 
but it is not preserved. In the case in which planets would 
be present in DG Tau, the most likely scenario would be the 
temporary Super-Earth capture into 3:2 commensurability 
by the Jupiter-like planet and then scattering. This case has 
been described shortly in Section [4. 2| Also for the RY Tau 
we do not expect the 2:1 commensurability. If in the inner 
disc of RY Tau there would be a pair of planets similar to 
that considered in our paper, then the Jupiter-like planet 
would migrate outwards and the Super-Earth inwards. The 
divergent migration would not induce any commensurabil- 
ity. 

We have also collected from the literature the observed 
properties of the circumstellar discs around other young 
stars. In Fig. [Tolwe have marked the location of several ob- 
served systems. Rox 44, UX Tau A and LkCa 15 have been 



observed and modelled by Espaillat et al. (20101. The data 



has been taken from their Table 3. The physical properties 
of Rox 44 and UX Tau A indicate the possibility of the oc- 
currence of the 2:1 commensurability. This is also the case 



of GM Aur studied by Hughes et al. (20091. Instead, the 



properties of these objects. Unfortunately, at radii smaller 
than 15 AU, which are the most relevant for our studies dis- 
cussed here, the surface density is model dependent and can 



disc in LkCa 15 is located well below the region where the 
resonant structure can be expected. If planets were present 
in this disc, the gas giant would migrate slower than the 
Super Earth and, since the Jupiter is on the external orbit, 
it would not be able to catch the Super-Earth and form the 
mean-motion resonance. The open circles are the locations 
of ten other discs studied by [Andrews et al.| ( |2010[ ) , namely 
Elias 24, GSS 39, DoAr 25, WAOph 6, VSSG 1, SR 4, SR 13, 
WSB 52, WL 18 and SR 24 S. They belong to the 1 Myr- 
old Ophiuchus star forming region. Two of those, namely 
DoAr 25 and SR 13, share almost the same properties and 
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Figure 11. An example of the 1:1 commensurability, which ends with the ejection of the Super-Earth from the system. The two upper 
panels show the semi-major axis evolution for both planets. The upper right panel is a zoomed out version of the short time interval 
presented in the upper left panel. The lower left panel presents how the Super-Earth eccentricity changes in time and lower right panel 
illustrates the behaviour of the resonant angle. 



are located very close to each other in our picture of the 
parameter space. 

Of course, the predictions coining from the observations 
are still somewhat vague due to the huge uncertainties in the 
determination of the disc structure parameters. The purpose 



of this Section was however to stress the fact that such pre- 
dictions are already possible and can be tested in the future 
with the advances of the observational techniques. 
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7 DISCUSSION AND CONCLUSIONS 

In this paper we have constrained the properties of proto- 
planetary discs which favour the locking of a Super-Earth 
and a Jupiter-hke planet orbiting around a Sun-like star in 
the 2:1 commensurability. The Super-Earth and the Jupiter 
are located on the internal and external orbits respectively. 
The results are illustrated in Figs. [T] and |10| In both fig- 
ures the medium grey region is where we expect the 2:1 
resonance to occur. In the low surface density regime (less 
or equal to 4 • lO"'') the final outcome is relatively easy to 
determine. The migration of both planets is rather slow in 
this region. We expect the 2:1 resonance for viscosities not 
lower than 4- 10~^. Below this threshold the Jupiter migrates 
slower than the Super-Earth and no commensurability can 
be achieved. On the other side, in order to get the 2:1 res- 
onance, the viscosity cannot be too high, because then the 
Jupiter will migrate outwards. This is the case in which the 
viscosity combined with the surface density through Eq. (pi 
gives accretion rates higher than 10~^ Mo/yr. In the high 
surface density regime (higher than 6 • 10~^) the situation 
is more complicated and the 2:1 resonance is allowed only 
in a narrow strip where the viscosity is varying in the range 
from 4 • 10~® till 10~^. In this regime the equilibrium gap 
profile is not empty. This is particularly true for relatively 
high viscosities. For that reason, putting the gas giant in 
the disc results in its very rapid migration. To avoid the 
artefact caused by the process of gap opening, we kept the 
Jupiter-like planet fixed, not allowing for its migration till 
the gap has reached the quasi-stationary state and stops 
changing significantly. In this way we have obtained the 2:1 
resonance for viscosities ranging from 4- 10~^ till 10~^. The 
results coming from numerical simulations have been com- 
pared with the analytical predictions summarized in Fig. [9] 
and with the current observations, see Fig. [To] 

The orbital migration of planets embedded in proto- 
planetary discs is a natural mechanism in which resonant 
configurations might occur. However, it is not the only way 
to form orbital commensurabilities. It has been suggested 
that the planet-planet scattering might be an efficient mech- 



anism leading to the resonant capture of planets ( Raymond 



et al. 20081. Our computational set up allows us to inves- 



tigate also this scenario in the situation when the gaseous 
disc is still present in the system. We have employed a simple 
procedure to locate both planets in a disc with high surface 
density (higher than 6 • 10~*) and let them migrate with- 
out waiting for the Jupiter to open a gap. Here we discuss 
in details one example of resonance obtained in this way. 
Our probe, which consists of a Jupiter-like planet and a 
Super-Earth, is located in a disc with surface density equal 
to 9 • 10"'* and viscosity value taken to be 8T0~^. As ex- 
pected, the Jupiter migrates extremely fast (with type III 
migration, |Masset fc Papaloizou (2003)), traveling one third 
of the distance from its initial location to the star in just 
100 orbits. Its orbit crosses the orbit of the slowly migrat- 
ing Super-Earth and we witness a resonant capture into 1:1 



mean motion resonance. This is shown in Fig. 11 (two upper 
panels). After the occurrence of the locking the Super-Earth 
stays very close to the Jupiter on the horseshoe orbit with 
an eccentricity of about 0.03 (Fig. 11 left lower panel). The 
resonant angle is shown in Fig. |11| (right lower panel). The 
1:1 commensurability created in the way described above 



lasted for roughly 700 orbits and after that the Super-Earth 
has been ejected from the system. 

Concluding, we have shown how the mean-motion res- 
onances observed in the already formed planetary systems 
can help in getting additional insight into the conditions oc- 
curring in protoplanetary discs in the early phases of the 
planetary system evolution. In other words, we have made 
an attempt to constrain the resonant properties of the plan- 
etary systems which may occur in observed protoplanetary 
discs. Assuming that the resonances are due to early stages 
of the evolution, we can answer the question on what should 
be the disc properties in order to attain the 2:1 resonance. 
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